arXiv:hep-ph/0503050v2 2 Aug 2005 


KEK-TH-1008 

KOBE-FHD-05-01 

HUPD-0501 


The Decay of Tau Leptons 
Produced in Neutrino-Nucleon Scatterings 

Mayumi Aoki^*, Kaom Hagiwara^’^ , 
Kentarou Mawatari^"^, and Hiroshi Yokoya"^’^^ 

^Theory Group, KEK, Tsukuba 305-0801, JAPAN 
‘^Department of Particle and Nuclear Physics, 

Graduate University for Advanced Studies, Tsukuba 305-0801, JAPAN 
^ Graduate School of Science and Technology, Kobe University, 
Nada, Kobe 657-8501, JAPAN 
'^Department of Physics, Hiroshima University, 

Higashi-Hiroshima 739-8526, JAPAN 
^Radiation Laboratory, RIKEN, Wako 351-0198, JAPAN 


Abstract 

Energy and angular distributions of the r decay products in the CERN-to-Gran 
Sasso n-r appearance experiments are studied for the decay modes t ^ -Kn and 
r ^ kvn {i = e or g). We find that the decay particle distributions in the laboratory 
frame are significantly affected by the r polarization. Rather strong azimuthal 
asymmetry of tt~ and i~ about the r momentum axis is predicted, which may have 
observable consequences even at small statistics experiments. 
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1 Introduction 

Neutrino oscillation physics is one of the most attractive held of current particle physics, 
and plenty of theoretical and experimental studies are revealing the amazing nature of 
the neutrino sector, such as their non-zero masses and the large mixings. We are now 
entering the stage of precise determination of the mass-squared differences including their 
signs and the mixing angles. The CP phase of the MNS (Maki-Nakagawa-Sakata) lepton- 
havor-mixing matrix pp also interests us. 

On the other hand, it is also important to hnd the direct evidence of the neutrino 
oscillation within the three generations. It can be achieved by detecting Ur appearance in 
the long baseline neutrino oscillation experiments with the initial beam. The CNGS 
(CERN Neutrino to Gran Sasso) long baseline experiments [21 with IGARUS [3] and 
OPERA |2] detectors aim to establish appearance by measuring the r lepton pro¬ 
duction events caused by the charged current (GG) interactions. They are now under 
construction and plan to start taking data from the year 2006. It is also expected that 
the GNGS experiments will improve signihcantly the current upper limit on ^ 13 , which 
is the smallest of the three mixing angles in the three-flavor MNS matrix, by measuring 

—*■ transition [Sj. The Ur GG events followed by the r —>• e decays contribute as 
background events to the signals of the > z/g events. Thus, the detailed analysis of the 
Ur GG events is important in the GNGS experiments. 

Since the r lepton has the large mass, m,- = 1.78 GeV, it immediately decays into 
several particles always including a neutrino (z^r)- For that reason, r production will be 
detected through its decay particle distributions. On the other hand, the decay particle 
distributions from r depend critically on its spin polarization. It is therefore important 
to consider the spin polarization of r in addition to its production cross section. 

Detailed discussions on the spin polarization of r produced in neutrino-nucleon scat¬ 
tering can be found in the recent paper IHj. There, the quasi-elastic scattering (QE), 
the resonance production (RES) and the deep inelastic scattering (DIS) processes were 
considered for the r production, and it was shown that the produced r^’s have high de¬ 
gree of polarization, and their spin direction depends non-trivially on the energy and the 
scattering angle of in the laboratory frame. 

In this article, we study the decay distributions from t~ leptons produced via the GG 
interactions, especially for the GNGS experiments. We consider the t~ production in the 
neutrino-nucleon scattering and its subsequent decays, for the decay modes t~ 7i~Ur 
and T~ —>• i~ueUr {i = e, /i): 



( 1 ) 


2 Tau production 

Let us start with the brief summary of the cross section and spin polarization of r 
produced via neutrino. One can hnd more details in Ref. [H]. 
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2.1 Kinematics and the formalism 

We consider t~ production by the charged current (CC) reactions off a nucleon target: 

Ur{k) + N{p)^T-{k’)+X{p’). (2) 

For the hadronic hnal states X, we consider three subprocesses; the quasi-elastic scattering 
(QE), the A resonance production (RES) and the deep inelastic scattering (DIS) processes. 
The four-momenta are parametrized in the laboratory frame as 

k = {Ey, 0, 0, Ey), 

p=(M, 0, 0, 0), (3) 

k' = {Er, Pr sin 9r, 0, Pr cos 9r), 

and the following Lorentz invariant variables are introduced 

Q^ = -q^ = -{k-k')\ (4) 

W‘^ = {p + qY. (5) 


Each subprocess is distinguished by the hadronic invariant mass W\ W = M for QE, 
M + < W < VFcut for RES. Went is an artificial boundary, and we regard that DIS 

process occurs in the regions of IF > IFcut^- We take IFcut = 1-4 GeV in this report. 

The differential cross section and the spin polarization vector of produced t~ are 
obtained in the laboratory frame as jB] 
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( 6 ) 
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(7a) 

(7b) 


(7c) 


where Gp is the Fermi constant and n = M^/{Q'^ + M^). s = {sx, Sy, s^) is dehned in the 
r rest frame in which the 2 ;-axis is taken along its momentum direction and the y-axis is 
along k X k', the normal of the scattering plane, in the laboratory frame. It is normalized 

^More detailed studies on r productions via the CC reactions have recently been reported in Ref. |7]. 


2 







as 1 ^ = 1/2 for pure spin eigenstates. are structure functions defined with the 

generic decomposition of the hadronic tensor, 

q) = -gf,uWi{p ■ q, Q^) + W 2 {p ■ g, Q^) - W^ip ■ q, Q^) 

. Tir / ^ 2 \ I P/j^qi' “i" qfiPu ttt / ^ 2 \ / q \ 

+ j^W^{p-q, Q ) + 2M^ W5{p-q,Q), ( 8 ) 

where the totally anti-symmetric tensor is dehned as 69123 = 1- These functions can 
be estimated for each process, QE, RES, and DIS, as follows. 


2.2 Hadronic tensor 


2.2.1 Quasi-elastic scattering (QE) 


The hadronic tensor for the QE process, Ur + n ^ t -|-p, is written by using the hadronic 
weak transition current as follows jHj: 

spins 


where 9c is the Cabibbo angle. is dehned as 

= {p{p)\J^^\nip)) = Up{p) T^Unip), ( 10 ) 

where T^ is written in terms of the six weak form factors of the nucleon, E// 3 , Fa, F^ 
and Fp, as 


P _ pV I '^^^J.aq i pV I ^pV I 
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{p + p1 
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+ —Fp 

M ^ 


75- 


( 11 ) 


We can drop F^ and F^ because of the time reversal invariance and the isospin symmetry. 
Moreover, the vector form factors F^ and F 2 are related to the electromagnetic form 
factors of nucleons under the conserved vector current (CVC) hypothesis: 


where 
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( 12 ) 

(13) 


with a vector mass My = 0.84 GeV and ^ = 3.706. For the axial-vector form factor Fa 
and the pseudoscalar form factor Fp, we use 


FA{q^) 


Fa{0) 

{1-qyMW 


FM") 


2 M 2 
ml — g 2 


FA{q^), 


(14) 


with Fa( 0) = —1.27 [Oj and an axial-vector mass Ma = 1.026 GeV ng. Notice that the 
pseudoscalar form factor Fp plays an important role for the polarization of r produced 
by neutrino because its contribution is proportional to the lepton mass and it has the 
spin-hip nature, although it is not known well experimentally; see Ref. m for details. 
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2.2.2 Resonance production (RES) 


The hadronic tensor for the A resonance production (RES) process, z/,- + n (p) —> r + 
A"*" (A++), is calculated in terms of the N-A weak transition current as follows [Hl IT^ 

ESI: 


lyRES 


4 ^ W‘^ - Ml + iWViW) 

c;mn« ^ ' ' 


(15) 


with the running width 


r(iE) = r(MA) 




(16) 


r(MA) = 0.12 GeV and A(a, b, c) = — 2{ab + bc + ca). The current for the 

process + n —>■ t~ + A~^ is parametrized as 


= {^'^{p)\JMip)) = '0“(pOr;,«Mn(p), (17) 


where ip°‘ is the spin-3/2 particle wave function and the vertex T^^ is expressed in terms 
of the eight weak form factors 456 ^® 


r = 

fia 


9^lc4-lpiQar.V , 9f^aP' ■ q-p'^qa^y , g^,aP■ q - p^^qa^v , q^q^r^v 

M + hP + AP 

, 9pa^-7pqa^A , 9paP' ■ q - p'^qa 

“T , ^ Oq -j- _ O4 -j- WaaOc -|- _ _r, Oe . 
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(18) 


By using the isospin invariance and the Wigner-Eckart theorem, we obtain another N-A 
weak transition current as (A++| J^|p) = v^(A+| J^ln). From the CVC hypothesis, Cq = 
0 and the other vector form factors 4 5 are related to the electromagnetic form factors. 
Assuming the magnetic dipole dominance [Hj, we have = 0 and CX = ■ For 

we adopt the modified dipole parameterizations 





(19) 


with CX (0) = 2.05 and My = 0.735 GeV. For the axial form factors, we use [TH] 


CX{q^) 


cm 



cm) 


m2 

m2 — g2 


cm). 


( 20 ) 


with ^5^(0) = 1.2 and Ma = 1.0 GeV. For and = 0 and = -\Ci 

good agreements with the data As in the case of Fp{q^) in the QE process, the 

pseudoscalar form factor G/^(g^) has signihcant effects on the r production cross section 
and the r polarization EH. 
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2.2.3 Deep inelastic scattering (DIS) 

In the DIS region, the hadronic tensor is estimated by using the quark-parton model; 

»r™(p.9) = E /f(21) 

q,q '' ^ 

where Pq = is the four-momentum of the scattering quark, ^ is its momentum fraction, 
and fq^q are parton distribution functions (PDFs) inside a nucleon. The quark tensor 
is 


q) =6{2pq-q-Q^ - ml,) 

X ‘^[-9llu{Pq ■ q) + 2,Pq^Pqu + Pq^q^ + q^^Pqjy =F ie/.uaf^Pq q^] ■ (22) 

The upper sign should be taken for quarks and the lower sign for antiquarks. We retain 
the hnal quark mass for the charm quark as rUc = 1.25 GeV, but otherwise we set 
rUg/ = 0. In the calculation, we used MRST2002 [T7] for the PDFs^. 

By neglecting both the nucleon mass and the initial quark masses, the following rela¬ 
tions are obtained: 

m2 

Wi{p-q,Q^) = F^{^,Q^), Wi=2,...Ap-q^Q^) = -(23) 

p-q 


Here, 


Fl = Y,fqM^Q^)^ F2 = 2j2Uq,qi^,Q"), 

q,q q,q 

Fs = 2j2 fqi^, Q")-2J2 ^4 = 0, F, = 2J2 fqM^ Q'), (24) 


<?><? 


where S, = Q'^/2p ■ q = x for massless hnal quarks [mq/ = 0), and ^ = x/X with A = 
/{Q^ + ^l') for q' = c. In fact, the differential cross section, Eq. ©, does not satisfy 
the positivity condition near the threshold with this naive replacement. We modify the 
VFi structure function as IFi = (1 -|- ^^)Pi in order to preserve the positivity [H]. 


2.3 Polarization of produced r at a fixed neutrino energy 

We summarize the cross section and spin polarization of t~ produced in neutrino-nucleon 
scattering for isoscalar targets at the incident neutrino energy = 10 GeV on the 
Pr cos Or-pr siu Or plane, where Pr and Or are the r momentum and the scattering angle in 
the laboratory frame. 

In Fig. m the differential cross sections da/dpz dpr, obtained from Eq. (jHl), are shown 
as a contour map^, where pz = Pr cos Or and pr = Pt sin Or. Only the contours of the DIS 
cross section are plotted to avoid too much complexity. Each contour gives the value of 

^We adopt a naive extrapolation of the parton model calculation for low in the W > IFcut region, 
by freezing the PDFs when < Qq (=1.25 GeV^). 

^It must be noted that this contour map differs from Fig. 11 of Ref. [H|, where da/dEr dcosOr is 
plotted. 
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Fig. 1: The contour map of the DIS cross section on the p,-cos sin 6^,- plane for 
the Vt-N —»• T~X process at Ey = 10 GeV in the laboratory frame. The kinematical 
boundary is shown by the thick curve. The QE process contributes along the boundary, 
and the RES process contributes just inside of the boundary. The r polarization are 
shown by the arrows. The length of the arrows give the degree of polarization, and 
the direction of arrows give that of the r spin in the r rest frame. The size of the 
100% polarization (P = 1) arrow is shown as a reference. The arrows are shown along 
the laboratory scattering angles, 9^ = 0°, 2.5°, 5°, 7.5°, and 10°, as well as along the 
kinematical boundary. 

the differential cross section in the unit of fb/GeV^; the outermost line is 1 fb/GeV^ and 
the innermost line is for 4 fb/GeV^. The kinematical boundary is shown by the thick 
curve. The QE process contributes along the boundary, and the RES process contributes 
just inside of the boundary. The contour map shows that the contributions in the forward 
scattering angles in the larger side are important. In that region, the cross sections of 
QE and RES are also large and comparable to that of DIS. 

The polarization vector s, Eq. dZD, of r is also shown in Fig. ^ The length of each 
arrow gives the degree of polarization (0<P = 2|s]<l)at each kinematical point and 
its orientation gives the spin direction in the r rest frame. The produced t~ have high 
degree of polarization, but their spin directions significantly deviate from the massless 
limit predictions, where all t~ should be purely left-handed. Since Sx of Eq. m turns 
out to be always negative, the spin vector points to the direction of the initial neutrino 
momentum axis. Qualitative feature of the results can be understood by considering the 
helicity amplitudes in the center of mass (GM) frame of the scattering particles and the 
effects of Lorentz boost from the GM frame to the laboratory frame. For instance, when 
we consider the neutrino-quark scattering in their GM frame, produced t~ has fully left- 
handed polarization at all scattering angle due to the V—A interactions and the angular 
momentum conservation. After the Lorentz boost, higher energy r’s in the laboratory 
frame preserve the left-handed polarization since those r’s have forward scattering angles 
also in the GM frame. On the other hand, lower energy r’s in the laboratory frame tend 
to have right-handed polarization since they are produced in the backward direction in 
the GM frame. Also note that < 0 holds always because the polarization vector points 
to the collision point in the GM frame. See more details in Ref. [Hj. Let us stress that 
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these features of the polarization of r play an important role in the following analysis. 

2.4 Polarization of produced t~ in the CNGS experiments 



£ (GeV) 

Fig. 2: The incoming Vj. flux with the three-neutrino model in the CNGS experiments. 
The neutrino oscillation probabilities are calculated for a set of parameters: = 

8.2 X 10“^ eV^ and 5m\.^ = 2.5 x 10“^ eV^, sin^ 26 ^ 12 , 23 , is = 0.8, 1, 0, respectively, the CP 
phase ^mns = 0°, and the matter density p = 3 g/cm^. 

In the CNGS experiments, beam is produced at CERN-PS, which is expected to deliver 
4.5 X 10^® protons on target (p.o.t.) per year. The beam is optimized for Ur appearance 
with a mean neutrino energy of about 17 GeV. Fig. |21 shows the expected Vr flux 

4>uAE.) = Y. Pu.^uAE.) (25) 

£=e, fj. 

at Gran Sasso with the baseline length of L = 732 km from GERN. Here (j)u^{Ey) are 
the initial vi fluxes {i = e, p) [2j and Py^^y^^Ey) are the vi —> Vr oscillation probabilities 
in the three-neutrino model. The fraction in the initial fluxes is less than 1%. 

The neutrino oscillation probabilities are calculated for a set of the three neutrino model 
parameters: two mass-squared differences 5m\2 13, three mixing angles 6*12,23,13 and the 
CP phase (5 mns in the MNS matrix [T], and the matter density p. The GHOOZ [THI and 
Palo Verde PI reactor experiments give the upper bounds on sin^ 26 * 13 , as sin^ 26*13 < 0.16 
for = 2.5 x 10“^ eV^. The values of 5m\2 and sin^ 26*12 are constrained by the 

observations of the solar neutrinos PI and the KamLAND experiment PI. and these of 
5m\^ and sin^ 26*23 by the atmospheric neutrino observation at Super-Kamiokande [22] and 
the K2K experiment [22]. No constraint on the CP phase has been given by the present 
neutrino experiments. In our analysis, we take the following values for the neutrino 
oscillation parameters: 

5m\2 = 8.2 X 10“^eV^, = 2.5 x 10“^eV^, 

sin^ 2012 = 0.8, sin2 2023 = l, sin2 20i3 = O, ^mns = 0°, (26) 

with the constant matter density of p = 3 g/cm^. Here we assume the so-called normal 
hierarchy. Because of setting sin^ 20 i 3 = 0, the neutrino oscillation probabilities are 
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Fig. 3: The contour map of the number of t~ production events on the cos 6^-p^-sin 6 t 
plane in the CNGS experiments. The t~ polarization are shown by the arrows. The 
length of the arrows give the degree of polarization, and the direction of arrows give that 
of the T~ spin in the t~ rest frame. The size of the 100% polarization (P = 1) arrow 
is shown as a reference. The arrows are shown along the laboratory scattering angles, 
Ot- = 2.5°, 5°, 7.5°, and 10°. The right hgure is an enlargement of the left hgure. 

approximately described by those in the two neutrino (z/^ and z/,-) model. See e.g.. Ref. |211 
for details. 

Taking into account the CNGS neutrino flux shown in Fig.|2l we show the distributions 
of events and polarization vectors of t~ on the Pr cos Or-Pr sin 6^ plane in Fig. IHl The right 
hgure is an enlargement of the left hgure to show the polarization vectors in detail for 
the important region of large cross section. The initial neutrino energy is integrated out 
with the incoming z/,- hux, of Eq. (|^ . whereas it is hxed at 10 GeV in Fig. [H 

The number of t~ production events for all the QE, RES and DIS processes are included 
in the contour map, where we assume 5 years with 4.5 x 10^® p.o.t./year of the primary 
proton beam and the 1.65 kton size detector, which are the current plan of the OPERA 
experiment |3]. Each contour gives a number of events per GeV®; the outermost line 
corresponds to 1 event/GeV®, and the innermost line is for 7 events/GeV®. The contour 
map shows that there are many events around Pt-= 10 GeV, and around 6r = 5°. As for 
the polarization vectors, the dependence on the energy and the scattering angle of t~ is 
rather smooth as compared to that in Fig. ^ because of the integration of the incident 
neutrino energy. As the right hgure shows, however, the direction of the r polarization is 
still non-trivial in the region which has many events. 


3 Tau decay 

3.1 Tau decay in the r rest frame 

Before turning to the r decay in the GNGS experiments, we give the formulas of the 
energy and angular distributions of the decay particles, especially 7t~ and i~, from the 












polarized t~ lepton in the r rest frame. The decay distribntions in the laboratory frame 
can be easily obtained by simple Lorentz transformation. Here again we note that the 
2 ;-axis of the rest frame, in which we calculate the spin polarization vector of r, is taken 
along its momentum direction in the laboratory frame. 

The decay distribution of tt~ via the decay mode r“ —>• vr”!/,- is given as 


1 dTjr 
Tr 






(27) 


where T.^ is the total decay width of r and = B{t ^ ttz /) = 11.06% [ 5 ] is the decay 
branching fraction. All the frame dependent variables with hat (") symbols are those in 
the r rest frame. The energy of vr is hxed as E.,t = ijni + ’>nl)/2mr, because of the 2-body 
decay kinematics. The angular distribution of 7i~ is dictated by the t~ polarization, and 
7T~ prefers to be emitted along the r polarization direction. 

Similarly, the decay distribution of i~ via the decay mode t~ —*• i~Pih'r = e,/i) is 

given as 


1 dVi 
Tr dEidCl. 
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1 2 


X PiEi 


3E„ 


2Et 


mj 2s-pe pe 
Ei Pe Ei 


E^ 
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(28) 


where .E^ax = (^r + Tnj)/2mr and the branching fraction of r —> iuu, B^, is 17.84% for 
i = e and 17.37% for £ = p, [Oj. The normalization function is /(p) = (1 — 8?/ -|- Sy^ — y^ — 
12y‘^\ny)/[1 + y^. The distribution of depends on its energy, and has a peak in the 
high energy. Furthermore, the high energy tends to be emitted against the direction 
of the T spin, although the impact is smaller than that of the tt case. It is meaningful 
to notice that, therefore, ti~ and tend to have the opposite preference of the r spin 
dependence on their angular distribution. 


3.2 Tau decay in the CNGS experiments 

In this section, we present our results of the decay particle distributions from t~ leptons 
produced by the CC interactions for the CNGS experiments. Main feature of our analysis 
is to deal with the proper spin polarization of t~ which is calculated for each production 
phase space, shown in Fig.|21 In order to show the effects of the r polarization on the decay 
distributions, we compare the results with unpolarized r decays and also with completely 
left-handed r decays. 

The events of the decay distributions for i = tt, £ are given by 


dN^ 

dEidVt\ 


rl 

=A / dEy (()y^{Ey) / d cos Or 

J Jc ■ 
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/ dEr 
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"" Tr dE-m, ’ 


(29) 
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Fig. 4: Energy distributions of vr” (thick lines) and e~, fi~) (medium-thick lines) in 
the decay of r“ produced in the neutrino-nucleon CC interactions for the CNGS experi¬ 
ments. 5 years running with 4.5 x 10^® p.o.t./year of the primary proton beam and 1.65 
kton detector are assumed. Solid, dashed, and dotted lines shows the energy distributions 
with the predicted r polarization, unpolarized, and purely left-handed cases, respectively. 
The estimated number of r“ production is also shown by a thin solid line with respect to 
the r energy (Er). 

where A is the number of active targets, is the maximum value of neutrino energy 
in the flux, = m,- -|-m^/2M is the threshold energy to produce r lepton off a nucleon, 
and the other integral ranges are given by 

Cmin{E,) = Vl + M/E, + M‘^IEI - mllAEl - 
E± {E,y, Or) = (& ± y/b'^ — ac) /a, 

where a = (Ej, -|- M)^ — E'^cos^Or, b = {Ey -|- M){MEy + ml/2), c = mlE/cos^ Or + 
{MEy + m^/2)^. Here, all the frame dependent variables are defined in the laboratory 
frame. Er and Or dependence of dT\/dE\(El\ appear by the Lorentz transformation from 
the T rest frame to the laboratory frame. Although we retain only the vr or £ momenta 
in the above formula, more exclusive measurements may be possible. For example, a 
well-developed emulsion detector used in OPERA experiment jTj can detect the kink of r 
decay and can measure the scattering angle of r with the accuracy of mrad, and also the 
r flight length is measurable. Therefore it would be possible to do more detailed analysis 
of r events, if there were enough statistics. 

Fig. m shows the energy distributions of vr (thick lines) and i (= e , p ) (medium- 
thick lines) decayed from t~ produced in the neutrino-nucleon CC interactions. We 
assume the same conhguration of the experimental setup as Fig. 01 i.e., 5 years running 
with 4.5 X 10^® p.o.t. per year of the primary proton beam and 1.65 kton size detector 
for the OPERA experiment |3]. (For the ICARUS experiment, 2.35 kton of liquid Argon 
detector mass and 10 years running are planned 01, so that more statistics are expected.) 
For each decay mode, solid lines show the distributions from the decay of r“ with the 
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Fig. 5: Azimuthal angle distribution of vr” (thick lines) and £~ (thin lines). The setup 
conhgurations, such as neutrino flux, detector size etc., are the same as Fig. HI Solid lines 
show the distributions from t~ with the predicted r polarization, and dashed lines show 
those with the unpolarized r case. The results of purely left-handed r are the same as 
those for the unpolarized r. 


predicted r polarization. For comparison, dashed and dotted lines show those of unpolar¬ 
ized T and purely left-handed r, respectively. The estimated number of the t~ production 
with respect to the r lepton energy, Er, is also plotted as a thin solid line. The results 
are calculated by using Eq. (Ei with 100% particle detection efficiency for simplicity. For 
the above parameters, 113 events of t~ are produced and 13 (20) of those decay into 7i~ 
(£~) mode. 

The 7r~ and £~ distributions have peak in the low energy region, and in this region 
the polarization dependence becomes large. As we pointed out in the previous section, 
the polarization dependence is opposite between vr” and £~, and is more significant in 
7r~ than in £~. In the peak region, the polarization dependence affect the distribution 
around 30% (15%) for the 7r“ (£~) decay mode. Expected statistics is rather small in 
the current design of the CNGS experiments. However, the likelihood probability of each 
event will be affected signihcantly by the r polarization effects. The characteristic feature 
of our prediction is that the produced t~ is almost fully polarized and that it has non-zero 
transverse component of the spin vector, namely of Eq. (ED- The observed patterns of 
the 7r~ and £~ energy distributions in the laboratory frame then follow from the energy 
angular distributions in the polarized r rest frame^. 

Fig. HI shows the azimuthal angle distribution of tt” (thick lines) and £~ (thin lines). 
The azimuthal angle (pi= 7 ry is given by dQi = d cos 6*i d(pi in Eq. ET?|l , and is measured from 
the scattering plane where (pi = tt/ 2 is along the p^, x direction in the laboratory frame, 
in which the z-axis is taken along the direction of the r momentum. Solid lines show 
the distributions from t~ with the predicted r polarization, and dashed lines show those 
from unpolarized r. The results of purely left-handed r are the same as those for the 

"‘We may note that we checked all our results by using TAUOLA (the Monte Carlo program to simulate 
decays of t leptons) ^H]. 
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unpolarized r. Since both unpolarized and purely left-handed r“ have zero component of 
perpendicular polarization, they give flat azimuthal distributions. The azimuthal angle 
distributions can be measured by tracking the trajectory of r leptons by emulsion detectors 
in the OPERA experiment, or by reconstructing the hadronic cascades from neutrino- 
nucleon scattering. As is the case of energy distribution, vr” and l~ decay mode show 
the opposite feature and polarization dependence is clearer on 7 r“ mode than t~ mode. 
At (p = 0 or | 7 r|, the dependence of the r polarization affects the distribution by about 
47% (16%) for the 7 r“ {i~) decay. Even though the number of event is limited, it may be 
possible to obtain a hint of such large asymmetries. 

Finally we comment on the contributions of the neutrino oscillation parameters to our 
results. The number of t~ production is very sensitive to the value of When we 

take Eq. P)|) but = 3 X 10 ^ eV^, about 50 more events of r are obtained. On the 
other hand, the produced r“ decreases (about 10 % maximum) for the larger (smaller) 
value of sin^26'i3 (sin^ 26 ^ 23 ). However those parameters do not change much the impacts 
of the T~ polarization on the energy and angular distributions of the r decay products. 


4 Conclusion 

In this article, we have studied the effects of the spin polarization of t~ produced in 
neutrino-nucleon scattering on the subsequent decay distributions. The calculation of 
the cross section and the spin polarization of t~ production processes, QE, RES and 
DIS, were reviewed and the decay distributions of t~ into 7i~ or £~{= e~, fi~) modes 
were considered. Taking into account the polarization of produced r“, we calculated the 
energy and azimuthal angle distributions of 7 r“ and i~ in the laboratory frame, for the 
experimental setup of the CNGS long baseline project, OPERA and ICARUS experiments. 

We found that the decay particle distributions in the laboratory frame are signihcantly 
affected by the t~ polarization. Rather strong azimuthal asymmetry of 7 r“ and i~ about 
the r momentum axis is predicted, which may have observable consequences even at small 
statistics experiments. 

Before closing, let us mention about the decay particle distribution from , although 
there is no plan to use beam in the CNGS experiments so far. In this case the strong 
azimuthal asymmetry predicted for the t~ decay is not expected because the transverse 
component of the polarization is rather small, as shown in Fig. 12 of Ref. jH]. 


Acknowledgments 

K.M. would like to thank T. Morii for encouragements and S. Aoki for valuable comments. 
H.Y. would like to thank J. Kodaira for encouragement and RIKEN BNL Research Center 
for helpful hospitality during his stay. The work of K.H. is partially supported by the 
core university exchange program of JSPS. 


References 

[ 1 ] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 


12 


[2] CNGS project home page, http://proJ-cngs.web.cern.ch/, 

[3] e.g. A. Rubbia, Nucl. Phys. B (Proc. Suppl.) 91 (2001) 223, 

see also ICARUS Collaboration home page, http://pcnometh4.cern.ch/, 

[4] OPERA Collaboration home page, http://operaweb.web.cern.ch/ 

[5] See e.g., V. Barger et al., Phys. Rev. D 65 (2002) 053016; P. Huber et ah, Phys. Rev. 
D 70 (2004) 073014. 

[6] K. Hagiwara, K. Mawatari, and H. Yokoya, Nucl. Phys. B 668 (2003) 364; Erratum- 
ibid. 701 (2004) 405. 

[7] K.S. Kuzmin, V.V. Lyubushkin, V.A. Naumov, Mod. Phys. Lett. A 19 (2004) 2815; 
K.S. Kuzmin, Seminar a KEK on Feb. 15, 2005. 

[8] C.H. Llewellyn Smith, Phys. Rep. 3 (1972) 261. 

[9] Particle Data Group (S. Eidelman et al.), Phys. Lett. B 592 (2004) 1. 

[10] V. Bernard, L. Elouadrhiri, and U. Meissner, J. Phys. G 28 (2002) Rl. 

[11] K. Hagiwara, K. Mawatari, and H. Yokoya, Phys. Lett. B 591 (2004) 113. 

[12] P.A. Schreiner and F. Von Hippel, Nucl. Phys. B 58 (1973) 333. 

[13] S.K. Singh, M.J. Vicente Vacas, and E. Oset, Phys. Lett. B 416 (1998) 23. 

[14] L. Alvarez-Ruso, S.K. Singh, and M.J. Vicente Vacas, Phys. Rev. C 57 (1998) 2693; 
and references therein. 

[15] E.A. Paschos, J.Y. Yu, and M. Sakuda, Phys. Rev. D 69 (2004) 014013. 

[16] M.G. Olsson, E.T. Osypowski, and E.H. Monsay, Phys. Rev. D 17 (1978) 2938. 

[17] A.D. Martin et al., Eur. Phys. J. G 23 (2002) 73. 

[18] The GHOOZ Gollaboration (M. Apollonio et al.), Eur. Phys. J. G 27 (2003) 331. 

[19] The Palo Verde Gollaboration (F. Boehm et ah), Phys. Rev. D 64 (2001) 112001. 

[20] B. Kayser, in Review of Particle Physics, Phys. Lett. B 592 (2004) 1, pp. 145-153, 
and references therein. 

[21] The KamLAND Gollaboration (T. Araki et al.), Phys. Rev. Lett. 94 (2005) 081801. 

[22] The Super-Kamiokande Gollaboration (Y. Ashie et ah), Phys. Rev. Lett. 93 (2004) 
101801. 

[23] The K2K collaboration (M.H. Aim et al.), Phys. Rev. Lett. 90 (2003) 041801. 

[24] M. Aoki et al., Phys. Rev. D 67 (2003) 093004; K. Hagiwara, Nucl. Phys. B (Proc. 
Suppl.) 137 (2004) 84. 

[25] S. Jadach, Z. Was, R. Decker, and J H. Kuhn, Gomput. Phys. Gommun. 76 (1993) 
361. 


13 



